REPORT No. 768 



THE FLOW OF A COMPRESSIBLE FLUID PAST A CURVED SURFACE 

By Carl Kaplan 



SUMMARY 

An iteration method is employed to obtain the flow of a com- 
pressible fluid past a curved surface. The first approximation, 
which leads to the Prandtl- Glauert rule, is based on the assump- 
tion that the flow differs but little from a pure translation. The 
iteration process then consists in improving this first approxi- 
mation in order that it will apply to a flow differing from pure 
translatory motion to a greater degree. The method Jails when 
the Mach number of the undisturbed stream reaches unity but 
permits a transition from subsonic to supersonic conditions 
without the appearance of a compression shock. The limiting 
value of the undisturbed stream Mach number, defined as that 
value at which potential flow no longer exists, is indicated by 
the apparent divergence of the power series representing the 
velocity of the fluid at the surface of the solid boundary. 

For small Mach numbers and for thin shapes, the results 
obtained by the iteration process agree with those obtained by 
the Poggi method. For higher values of the stream Mach 
number less than the critical value, numerical calculations are 
in agreement with the results obtained by von Kdrmdn by means 
of the hodograph method. For values of the stream Mach 
number higher than the critical value, the iteration process 
yields some information about the region of flow comprised 
between the critical stream Mach number and the limiting 
stream Mach number. 

INTRODUCTION 

When a body is held fixed in a compressible fluid moving 
at a uniform speed less than, but comparable with, that of 
sound, there may be a region near the surface where the 
velocity of the fluid relative to the body exceeds the local 
velocity of sound. The flow in such cases may be perfectly 
regular with no indication of shock waves. Several such 
types of flow have been described by Taylor (reference 1) 
and, more recently, by Gortler (reference 2). In connec- 
tion with this type of flow it is important to know when the 
irrotational motion ceases to be possible. It is certain that 
irrotational motion is no longer possible as soon as the 
Mach number of the undisturbed stream reaches a definite 
value, always less than unity, which depends on the shape 
of the body. In the past, compressibility shock has often 
been assumed to occur when the maximum velocity of the 
fluid at the surface of the body equals the local velocity of 
sound; however, the papers of Taylor and Gortler (refer- 
ences 1 and 2) question whether this is correct for the first 
appearance of a shock wave. In addition, a recent paper by 



von Karman (reference 3) suggests that the envelope of the 
Mach lines in the supersonic region of flow probably intro- 
duces the first shock wave in the flow. The stream Mach 
number at which the envelope of the Mach lines first ap- 
pears may thus be identified with a limiting value of the 
Mach number. 

The present paper treats the flow of a compressible fluid 
past a curved surface by means of an iteration process based 
on that of Ackeret (reference 4). The boundary was so 
chosen as to conform with the requirements of the method; 
namely, no stagnation points and small variation of the local 
velocity from that in the stream. The process, further- 
more, permits values of the stream Mach number ranging 
from zero to the neighborhood of unity. The method 
proves to be quite laborious when more than two stages in 
the iteration are demanded; but, because of the importance 
of the problem, it has been thought worth while to perform 
the third step. Most of the details of the calculations have 
been relegated to appendixes in order not to disturb the 
continuity of the main ideas. (The. equations in the appen- 
dixes have been assigned numbers prefixed by letters 
denoting the appendix; for example, equation (A-3) is the 
third equation in appendix A.) 

THE ITERATION PROCESS 

The fundamental differential equation governing the flow 
of a compressible fluid is 

f--> <*-<> |h- (&t&)-o 0) 

where 

X, Y rectangular Cartesian coordinates in plane of profile 
u, v fluid velocity components along X- and F-axes 
c local velocity of sound 

The condition for irrotational motion is that 

bu by 

bT~bX 

and leads to a velocity potential <t> defined by 

b<f> 



u= 



v= 



bX 



(2) 
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If the body is held fixed in a uniform stream of velocity U, 
the relation between the local velocity of sound c and the 
velocity of the fluid •yV+i' 2 is given for adiabatic processes by 



1- 



u 2 + 



(3) 



where 

Co velocity of sound in undisturbed stream 
7 ratio of specific heats at constant pressure and constant 
volume 

M Mach number of undisturbed stream (U/co) 

With the introduction of a characteristic length s as unit 
of length and the stream velocity U as unit of velocity, the 
various quantities thus far denned can be rendered non- 
dimensional. Thus X, Y, u, v, and <f> denote, respectively, 
the nondimensional quantities X/s, Y/s, u/U, v/U, and <t>/Us 
while c and Co retain their original meanings. By use of 
equation (2), equations (1) and (3) then become, respectively, 



and 



&4> 



-2M 



d<j> d<j> 



a»4 

bXbYbXb 




il^l-C^+c 2 )] 



(5) 



Let t denote a characteristic parameter of the shape, such 
as the thickness coefficient; then, the following expansions 
are assumed: 



(6) 



When these expressions for u and v are introduced into 
equation (4), together with the expression for c?/cq s given by 
equation (5), and when the coefficients of the various powers 
of t are equated to zero, the following differential equations 
for fa, fa, fa, . . . result: 



<7> 



n 7i,fA &fa &fa_ 

+ 2 ^aX aF bXbY + bY bXZY + bT bXZYJ f w 



These differential equations may be put into more familiar 
forms by introducing a new set of independent variables x 
and y, where 

x=X 1 

J (10) 
y=J[=WY) 

Thus, for M<1, equation (7) is transformed into a Laplace 
equation and equations (8) and (9) are transformed into 
Poisson equations. The solution of equation (7) yields the 
well-known Prandti-Glauert rule, whereas the solutions of 
equations (8) and (9) provide higher approximations to the 
flow of a compressible fluid and thus will apply for larger 
departures from an undisturbed uniform flow. 

The procedure to be followed in solving equations (7), 
(8), and (9) is very simple in principle. The first step is to 
obtain an expression for the velocity potential of the incom- 
pressible flow past the chosen boundary and to express it 
as a power series in the thickness coefficient t. Then the 
solution for the first approximation fa to the compressible 
flow is easily obtained by analogy from the coefficient of 
the first power of /. The second and third approximations 
fa and fa are obtained by solving equations (8) and (9). 
The boundary conditions — that the flow is tangential to the 
solid boundary and that the disturbance to the main flow 
vanishes at infinity — are satisfied to the same power of the 
thickness coefficient t which is involved in the expression for 
the velocity potential fa 

FLOW PAST A CURVED SURFACE 

The solid boundary chosen for use in this paper is a sym- 
metrical shape with cusps at both the leading and the trailing 
edges, thus insuring no stagnation points in a uniform flow 
parallel to the axis of symmetry. Appendix A contains the 
derivation of this shape and also the solution for the flow 
of an incompressible fluid past it. Appendixes B, C, and D 
contain the detailed calculations for fa, fa, and fa, respec- 
tively. The final expression for the velocity potential 4> 
takes the following form: 

<7>=cosh£ cos v+-£p (3e~ f cos 17— e -3 * cos Srj) 
+P(Fi cos ri+F 3 cos 317+^9 cos 5ri) 
+f cos 17+ 63 cos 3i7+<?s cos 5ri+6 7 cos 7v 



-e _7£ cos 



7*)]+ • • 



(11) 



where 

£, r] elliptic coordinates related to the Cartesian coordinates 
x and y by the equations 

z= cosh £ cos i) 
y= sinh £ sin i\ 



(12) 
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Fi, Ft, F B functions of £ and of M given by equations (C-18) 
Qi, Q s , O c , G 7 , 6 0 functions of £ and M given by equations 
(D-16), (D-17), (D-18), (D-19), and (D-22), respectively 

Equation (11) represents the solution of the fundamental 
differential equation (1) that satisfies the boundary condi- 
tions at the surface of the body and at infinity, insofar as 
the terms inclusive of the third power of the thickness 
coefficient t are concerned. The coefficients of the various 
powers of t are exact and are valid for all values of the Mach 
number M from zero up to but not including unity. On 
the other hand, the method of Poggi yields the components 
of the fluid velocity in the form of power series in M 2 , the 
coefficients of which are exact and valid for the entire range 
of values of the thickness coefficient t. Appendix E con- 
tains the solution of the problem of this paper by the method 
of Poggi as far as the M 2 terms are concerned. 

Because the iteration process and the Poggi method yield 
solutions of the same equation (1) in the form of power 
series in t and in M 3 , respectively, the two methods must 
agree in the range common to them; that is, the iteration 
expression for the fluid velocity at the solid boundary, 
expanded according to powers of M and with all terms con- 
taining powers of M higher than the second neglected, must 
agree with the corresponding Poggi result, expanded accord- 
ing to powers of t and with all terms containing powers of 
t higher than the third neglected. This calculation is shown 
in detail in appendix P. 

NUMERICAL APPLICATIONS 

Calculations are now made for the velocity distribution 
at the surface of a bump — that is, a member of the family of 
shapes derived in appendix A — for several values of the 
Mach number. Because terms involving powers of the 
thickness coefficient t higher than the third have been 
neglected throughout the present paper, the fluid velocity g 
should be expressed in the following form: 



2=l+aif+a s i ! +a s f+ . . . 

From equations (F-6) it follows easily that 

3 • 

<h—~ 2p cos 2a 



(13) 



cos 4c 



a*=-^ (l-« + (G r 1 +3G f ,+5Gf s +7G' 7 ) 0 



(l-/S)+2(3<? s +5<? 6 +7^)o] cos 2a 
f2(5<?H-7<? 7 ) 0 ]< 
14(&)o] . 



+ 



4+5/S+3/3 2 

33+82/3+14/S 2 
J8 



| cos 4a 
I cos 6a 



(14) 



The expressions for (Gi) 0 , (G s ) 0) (G 6 ) 0 , and (G 7 ) 0 are 
given at the end of appendix F, and table rV shows the 



calculated values for M=0.50, 0.75, 0.83, and 0.90. Table V 
gives the calculated values of Oj, a^, and a 3 at various 
positions along the profile for M=0.50, 0.75, 0.83, and 
0.90. With 2=0.10, the expressions for the maximum fluid 
velocity at the surface can be written as follows: 

t ■ C f 

2„„*=1+0.17321+0.02274+0.00344 for M=0.50 

2m«x=l+0.23026+0.05435+0.02188 for M=0.75 

g* a *=l+0.26893+0.11014+0.07240 for M=0.83 

tf««x=l+0.34412 +0.19268 +0.43784 for M=0.90 



(15) 



An examination of the foregoing series shows that q max will 
probably diverge for some value of M in the neighborhood 
of 0.83. This value of M marks the limit of irrotational 
potential flow and probably indicates the first appearance of 
a compression shock at the solid boundary. Farther on 
in this section a rule will be formulated for estimating this 
limiting value of the stream Mach number. 

The velocity distribution for a profile of thickness coeffi- 
cient 2=0.10 is calculated by means of table V and equa- 
tion (13). Table VI lists the values of g for Af=0.50, 0.75, 
and 0.83 and figure 1 shows the corresponding graphs. 
The broken curve represents the velocity distribution for 
(=0.10 and M=0.50 calculated according to the Poggi 
method. (See table ILT.) The curves of figure .1 show the 
agreement between the values of g calculated by moans 
of the Poggi and the iteration methods for Af=0.50 and 
also the gradual change in curvature of the velocity- 
distribution curves in the neighborhood of the leading and 
the trailing edges as the stream Mach number is increased. 

The critical Mach number, denned as that value of the 
stream Mach number at which the local fluid velocity 
first attains the local speed of sound, is calculated as follows: 



In equation (5), (w'+o 2 ) is put equal to jj& or 



2- S =^l(]|5+T-l) 



(16) 



Table VII lists values of g„ for various values of the Mach 
number M. From equations (14) with ct=^ and the expres- 
sions for ((?i) 0 , (#3)0, (6^)0, and ((? 7 ) 0 , it follows that the 
maximum velocity at the boundary is given by 



*3 4+3/3- 3 
/3V + 8 p 



"I 



2«ai=l+^ f+j~| (r+1) 
"f [~ m b+O ^ffl+r«0.6C-lOZ» 

+ (f+l) (-64+335-2.45<7+3.3l>)] [<•+... (17) 
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Fioubk L— Velocity distribution at the surface of a bump, f-CUO, for several values of the Mach number. 



Table II gives the values of A, B, G, and D obtained from 
equation (D-3) for various values of the Mach number. 
Table VII lists the values of g^ for f=0.10 calculated by- 
means of equation (17). Values of g,^ calculated by the 
Poggi method are also given in table VII. For low Mach 
numbers the approximate values for obtained by the 
two methods agree, but for high Mach numbers the Poggi 
method yields values that are too low. Figure 2 shows the 
graphs of q„ and of q nax as functions of the Mach number. 
The intersection of the two curves gives the critical Mach 
number. The iteration method (solid line) yields the value 
Mer= 0.742, whereas the Poggi method (broken line) yields 
the value M CT =0.788. 

The value of q a for M=0.83 obtained from equation (16) 
is 1.1731. The last column in table "V 1 1 1 shows the values 

of qJicr for the shape f=0.10 with M =0.83. Values of -^->1 

So- 
designate a region in which the flow is supersonic, whereas 



The 



values of ^"<Cl characterize the subsonic region. 

go- 
supersonic region of flow is symmetrical with regard to the 
F-axis and corresponds to a calculation given in reference 1 
for the flow through a nozzle, in which a similar supersonic 
region of flow was found at the wall at the narrowest cross 
section of the tube. 

In order to find the extent to which the supersonic region 
penetrates the flow, it is sufficient to utilize for the calcula- 
tion only the terms inclusive of C, since the series for g 
converge rapidly away from the solid boundary. Thus, 
along the F-axis, the following expression for g is obtained 
from equations (F-3) and (F-4) : 



-D(8€-«-5«- 3 f+3e- 5 «)]Jf s + . 



(18) 
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Fioube %— Marimnm velocity at the surface of a bump, <«=0J0. 



Then for <=0.10 and M=0.83, 

2=l+0.40316e- s f+0.02753e- 4 « 
+ cdhl (-°- 02762 «" t +°- 01242 «" 3f -0-0 1 520e- 6 «) 

The value of £ for which 2=2 CT =1.1731 is 0.38. The corre- 
sponding value of Y obtained from equation (B-16) is 

Y=\ sinh £ 
p 

=0.720 

■The supersonic region of flow thus extends into the fluid a 
distance equal to almost seven times the maximum height 
of the bump. By use of several chosen values of 17, the 
constant velocity profile 2=2o-= 1.1731 can be plotted by 
means of equations (F-3) and (F-4). In figure 1 the lower 
broken curve represents this profile. The region inside the 
■ profile is completely supersonic and therefore contains real 
Mach fines. The region outside the profile is everywhere 
subsonic and therefore the Mach lines are imaginary. 

The pressure coefficient O v ^i is obtained from the 
expression 



r/7-1 



(19) 



where 

\poTP 

and 2 is the velocity of the compressible fluid, referred to 
the velocity U of the undisturbed stream. 

Since, throughout this paper, terms involving powers of t 
higher than the third have been neglected, C pJ s should be 
expressed as a power series in t. Thus, if g from equation 
(13) is substituted into equation (19), it follows easily that 

0„, m= —2(ht+ [— (a 1 J +2a 2 ) +a 1 i M z \ i 3 

+[-2(a 1 a 1 +a3) + (a 1 s +2a 3 )a 1 il4 5 - 2 ^a 1 3 M < ]< 3 + . . . (20) 

With the help of table V, values for G PwU along the profile 
can be easily calculated for the case f=0.10 and _M=0.83. 
Table IX shows these values of G PM together with correspond- 
ing values calculated according to the Prandtl-Glauert 
rule and the von K ton an method. The Prandtl-Glauert 
rule is 

Cp.o 



C Pw M—- 



and the relation obtained by von Karman (reference 3) is 



Gp.M — 



Figure 3 shows the graphs of the various calculated results. 
It is observed that the results of the Prandtl-Glauert and 
the von Karman methods differ considerably from the 
results of the iteration method. The reasons for these differ- 
ences are that the Prandtl-Glauert approximation, though 
valid for Mach numbers in the neighborhood of unity, should 
be utilized only for very thin shapes ; whereas the von Karman 
method, though applicable to any reasonable shape, is no 
longer valid for Mach numbers beyond the critical value. 

As final numerical applications of the results of this paper, 
the maximum values of the negative pressure coefficient 
— (Cj>,if)moi, the critical pressure coefficient — (C p ^r) CT , and 
the limiting pressure coefficient — (Gpj^iim are calculated 
for various values of the thickness coefficient and of the 
stream Mach number. 

The maximum values of the negative pressure coefficient 
for various values of t and of M are obtained by means of 
equations (17) and (20). Table X shows values of a u da, 
and 03 defined, respectively, as the coefficients of f, f 2 , and f 3 
in equation (17). The corresponding values of G pM , ob- 
tained by means of equation (.20), are listed in table XI, 
together with values of O pJJ calculated by the von Karman 
method. Pigure 4 shows the variation of {G P ji)m<a with 
Mach number for several values of t. 

The critical pressure coefficient (Gp^a- is calculated by 
means of the following expression obtained by substituting 
for 2 s in equation (19) the expression for from equation 
(16): 
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-Pressure distribution at tha surface of a bump, <=0J0, for 3f=<UB. 



Table XII gives the values of {.Gvasla calculated by means 
of this equation, and figure 4 shows the corresponding graph. 
The intersections of the {Oj,^) max curves with the {C vM )„ 
curve yield the critical values of the stream Mach number 
which are listed in table XII. 

As noted once before, for a given value of t, the series for 
2 (equation (13)) apparently diverges for a definite value of 
the stream Mach number. It is reasonable to assume that 
this value of the Mach number marks the limit of irrotational 
potential flow and also probably indicates the first appearance 
of a compression shock at the solid boundary. Equation (17) 
for q maT can be used to estimate the limiting values of M 
according to the following rough criterion. By means of 
table X, expressions for q max , in the form of power series in t, 



-2.6 



-2A 



-2.0 



J 



-I. -4 
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Fiotoe 4.— Mairlmnm pressure coefficient as a function of the Mach number. 

critical pressure coefficient at which local fluid speed equals local speed of sound. 
(C„u)iim limiting value of pressure coefficient beyond whloh flow without shook 

cannot exist. 

i.C„u)tu absolute limiting pressure coefficient; that Is, corresponding to a vacuum. 



can be obtained for any value of the stream Mach number in 
the range 0^M<1. For a given value of the thickness 
coefficient t, a series for q max can then be found so that a term- 

°° 1 

for-term comparison with the harmonic series S - yields a 

decreasing sequence for the ratio of corresponding terms. 
The value of the Mach number thus obtained is chosen as the 
limiting value of the Mach number. Table XII lists both the 
values of (Opj^nm and the corresponding values of M for 
different values of t. Figure 4 shows the curve connecting the 
limiting values of O v ,u. The region between this curve and 
the (C Pt n)er curve represents the supersonic range without 
compression shocks. It is emphasized that the mathematical 
procedure outlined in this paragraph is a highly speculative 
one. 
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The maximum speed that can be attained by an adiabatic 
fluid is equal to the speed when p=p=c=0 and is given by 
Bernoulli's equation 

2 

2 a6 . 2 =l+ (7 _ 1)M s 
Values of the pressure coefficient 

{0 P , it) ab ,= (22) . 

are listed in table XTtT for various values of the undisturbed 
stream Mach number M, and figure 4 shows the correspond- 
ing curve. The region between this curve and Hi&(0^^)ur, 
curve represents the supersonic range with compression 
shock. On and beyond the {O v ,u) ab , curve the adiabatic fluid 
ceases to exist; that is, absolute vacuum prevails. 

In conclusion, it may be remarked that the results of 
von Karm&n, shown by the small circles in figure 4, are 
obtained independently of any assumption (MJncerning the 
shape of the solid boundary; whereas the results of this 
paper were obtained for a specific family of shapes. As 
shown by the curves of figure 4, nevertheless, the results of 
this paper agree with those obtained by von Karman's 



method. This agreement has some justification, for the 
values of (Op iir )max and M„ depend mainly on the dimensions 
of a body — that is, thickness coefficient — and not on its 
shape. The hodograph method as employed by von Karman, 
however, yields results that cease to be valid when the 
local Mach number equals or exceeds unity; whereas the 
iteration method utilized in the present paper is valid for 
local Mach numbers greater than unity and, for the family 
of shapes chosen, yields some information with regard to a 
supersonic region of flow without shock. The limiting 
value of M for a given shape, beyond which supersonic 
flow without shock does not exist, appears to depend on 
the convergence of the power series in t for the velocity. 
Although only a few terms of the series have been obtained 
and therefore the limiting value of M cannot be given 
precisely, nevertheless it is believed that a reasonable esti- 
mate of the value for M lim can be made by the comparison 
test with the harmonic series. 



Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langlby Field. Va., September IS, 194S. 



APPENDIX A 



THE INCOMPRESSIBLE FLOW PAST A CURVED SURFACE 



In the search for a shape which satisfies the conditions that 
it be thin and that it possess no stagnation points, the first 
thought is of a straight-line segment. It is well known that 
a straight-line segment of length 4c is obtained from a circle 
of radius c by means of a JoukowsM transformation. If Z 
denotes the plane of the segment and Z' the plane of the 
circle, then 

e 2 



Z=Z'+^, 



(A-l) 



The singular points of this transformation are determined by 
the equation 

which shows zeros at Z'=±c. In order to raise the top 
surface and lower the bottom surface of the line segment, it 
is necessary only to place two additional zeros at Z'=±id 
where (2<c. Analogous to equation (A-2), 

S-( i +#)( i -l)( i +D0-D < A - 3 > 



Then, on integration of equation (A-3), it follows that 



Z=Z' 



zz n 



(A-4) 



The parametric equations of the shape in the Z-plane corre- 
sponding to the circle of radius c in the Z'-plane are obtained 
by substituting Z'=ce K in equation (A— 1); thus, 



d? 

X=2c cos 6 — ^ (3 cos 6 — cos 30) 
Y~ (3 sin 5 - sin 39) 



(A-5) 



The family of shapes given by equations (A-5) includes, on 
one hand, the straight-line segment with d=0 and, on the 
other hand, a shape having four cusps symmetrically, placed 
with respect to the coordinate axes with d=c. For 0<d<c, 
the slope dY/dX is zero for 0=0, tt/2, and ■sr. The shape thus 

has cusps at X=±2c (l— Y=0, and the maximum 

4 d 1 

and minimum points are at X=0, Y=g — > and X=0, 
4 d 3 

Y=— g —> respectively. 

The complex potential for a circular cylinder of radius c, 
fixed in a stream of uniform velocity U in the positive 
direction of the real axis, is given by 



F=u(z'+£) 



(A-6) 



The complex velocity past the corresponding shape in the 
Z-plane is 

. dF dF dZ' 
u ~ w= dZ=dZ' dZ 

or 

, 2 . . dF dF dZ' dZ' 
?= u +* = dT?dZ>dT d 2 

By means of equations (A-4) and (A-6), it follows that 



0 



~^d 2 



1 - Z /2 %,<) 



~W Z 

and at the surface of the profile, where Z'=ce <s , 

IP 



1+2 ? cos 26+(' P ^ 



(A-7) 



It will be convenient to consider F, Z, and Z' as non- 
dimensional quantities. Thus, in the plane Z' the unit of 
length is the radius c of the circle; whereas, in the plane Z, 
the unit of length is the semichord s of the shape. Then 
Z', Z, and F denote Z'\c, Z/s, and FJUs, respectively. Also, 
the ratio (F/c 2 is designated by e. 

According to equations (A-5), the semichord s is given by 



and the thickness coefficient t by 

2e 



f= 



3-6 



With the introduction of these new designations, equations 
(A-4), (A-5), (A-6), and (A-7) become, respectively, 



4 



1—t 1 , t 1 
2 Z' + 4Z ,J 



X= cos 0 — ^ (cos 0 — cos 30) 



(3 sin d — sin 30) 



^(z< + i.,) 

^ = l+2e cos 20+e 2 



(A-8) 
(A-9) 

(A-10) 
(A-ll) 
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p 0 static pressure in undisturbed stream 
p density of fluid 
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As a numerical example of the use of equations (A-9), 
(A-10), and (A-ll), table I gives the coordinates for the 
shape £=0.10 together with the velocity and pressure dis- 
tributions along the profile. The pressure coefficient C v is 
calculated by means of Bernoulli's equation: 



where 

p local static pressure 



Figure 1 shows the curves of the shape and the velocity 
distribution, and figure 3 shows the graph of the pressure 
distribution. Because the body is placed symmetrically 
with respect to the undisturbed stream, the flow is identical 
with that over a solid boundary composed of the X-axis 
from X= oo to X=l, the upper surface of the shape with 
O^d^Tr, and the X-axis from X=— 1 to X= — co. This 
boundary is called a bump. 



APPENDIX B 

THE COMPRESSIBLE FLOW PAST A BUMP 



Before proceeding with the iteration process, developments 
for <f> and Y in positive, integral powers of the thickness 
parameter t must be obtained. From equation (A-8), 



f=4 



and from equation (A-10), 



(B-l) 



(B-2) 



By means of a Taylor expansion in the neighborhood of 
t=0, it is possible to express .Fas a power series in t in which 
the coefficients are functions of Z. Thus, according to equa- 
tion (B-l), t=0 when 



or 



Z'=z+yz*-i 

where the positive sign is taken with the radical because the 
points at infinity of the Z- and Z'-planes must correspond. 
Now 

r-W~+t(*g)„+ ■ • - 

If, on the right-hand sides of equations (B-l) and (B-2), Z is 
considered constant and Z' is considered the parameter, it 
follows that 

F=Z+[|z-(Z 2 -l)(Z-VZ 2 ^T)]f+ • • • 

Since F=<f>+i^, where ^ is the stream function, 

4,=X+l[x-(Z 1 -l){Z~ VZ^D-CZ 5 - 1)(Z- V2*=l)~] *+ 

21 J (B-3) 

In a similar manner, equations (A-9) con be written as 
follows: 

X— cos 0 



' sin 2 6 cos 6 
Y= — (X— cos (9) tanflj 



(B-4) 



If X is considered constant and 0 is considered the parameter, 
a Taylor expansion yields 

F-a-XT^-SX^-i-! (SX*-3X*)P- . . .] (B-5) 

DETERMINATION OF 

By means of the transformation x=X, y=f)Y, where 
j3=V 1— M*, equation (7) for fa becomes 



(B-6) 



By analogy with the coefficient of t in equation (B-3) for 4>, 
it is suggested that 



*,=*[*- (2 s - 1) (z- V?=T) - (?- 1) (5- V^l)] (B-7) 

where the coefficient k is determined by means of the bound- 
ary condition 

00 dY 
bXdX~W 

or 



(B-8) 



and where the boundary, obtained from equation (B-5), is 
given by 

y=P(l-x*)*' i [<-3xV+| (8x i -Zx r )f+. . .] (B-9) 

It is clear that the boundary condition need be satisfied only 
to the same degree in t as is involved in the development of 
the velocity potential <j>. Thus, 



4=x+k[s-&-l)(g--J2^)-&—i)(2—-jW=i)]t+ • • • 

(B-10) 



(B-ll) 



^=i(||-||)=-^(-V?=i)-2^-VP=i) 

•V^^(s-V? zr l) + V^ r: T(^-V^ =r l)K+ • • • (B-12) 
^=-3ftrVra 1+ . • - (B-13) 
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Suppose that 2= cos a. Then, insofar as the terms of the 
first power in t are concerned, it follows that, on the boundary, 

2= cos a+iflt sin 3 a 

2= cos a—ifit sin 3 a 



•vV— l=i sin a(l— iftt sdn 3 a) 
V i 2 — 1 = — i sin a(l +i/3f sin 3 a:) 
Then, on the boundary, 

|^=1— 3fa cos 2« 

|2=-3fosin2« 
by 

g» — lo^s. 

The boundary condition, equation (B-8), then yields 



Therefore 



K 2/S 



+^3 [s- (2*- 1) (2- V?=l) -(?-!)(£- V?~ 1 )]*+•• - 



This expression for 4> can be put into a simple form by means 
of the transformation 



where 
Also, 



z= cosh f 

i=cosh I cos 77] 
y=sinh £ sin 17 J 



(B-15) 



(B-16) 



Then from equation (B-14) it follows that 
<£o= cosh £ cos 17 ' 

♦ 

#1= jjg (3e _t cos ij— e~ 3£ cos 

INVERSION OP EQUATIONS (B-16) 



(B-17) 



The relationship between the rectangular Cartesian 
coordinates X, T and the elliptic coordinates £, 17 is obtainod 
as follows: 

Invert equations (B-16) and solve for £ and rj; thus 



and 



Solve for sinh 2 £, 



(cosh f) + (sinh {) ~ 1 

(-*-) 2 -(- y -) 2 =l 
\cosh £/ \sinh f / 



2sinh^=-6+V& 3 +4y 2 



and solve for sin* 77 



2 sm 2 ,=6+7^+4? 

where 

6=1 -W+tf) 
By means of the transformation 

x=X 
y=BY 

it follows that 

2sinh 2 £=-6+V6 4 +4/3 i r j 



where 



APPENDIX C 

DETERMINATION OF fa 



The differential equation (8) for <fo in terms of the 
coordinates x and y becomes 



by 2 



By means of the symbolic relations 
be b2 55 

y_y 0 b 2 d 2 

b* 2 d2 2-1 "- J b2dl" t "d? 



2sin 2 7 7 =6+V6 2 +4^r 2 
6=1-(Z 2 +/3 2 F ! ) 

by \b2 52/ 



b 2 



b 2 0 s 



by 2 bz^^bi bl 5 



b 2 ■ / b 2 b 2 \ 
bxby -1 ^ bSV 



and the transformation 



2= cosh f 



(B-18) 



equation (C-l)_can be expressed in terms of the complex 
variables f and f as follows: 
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sinh f siuh J 5f /3 2 



»4*__i-PU, +1) , nfl -/ 1 ^i. i SM/ i &V, coshr a^i , i a»fr coshjd&N 



1 foi 1 / 1 d 2 ^ cosh f dfo 1 d s ^ , coshj d<fr 



Now, equation (B-17) for <fo can be rewritten as follows: 



* 1= £j3 (3e _r -e- 3r +3e-f-e- 3 0 



(C-3) 



The substitution of this expression for <£i into the right-hand 
side of equation (C-2) then yields the following differential 
equation for cfe: 

W + W^ ^U(7+D + (3-T)^](e^- e -3 £ ) cos, 

+ (7- 1) (1 - /S 2 ) (c- { f cos 3q-e- 3 « cos 5ij) } (C-4). 

The right-hand side of this equation suggests a solution of 
the form 

<h=Fi(0 cos n+F,(Q cos 317+^(5) cos 5 V (C-5) 

If this expression for <fo is inserted into equation (C-4) and 
the coefficients of cos 17, cos 3i7, and cos 5?? are equated on both 
sides of the equation, the following differential equations for 
Fx ft), and F&) are obtained: 

W --^fH^rWl) + (3-7) m (e-*-e-'i) (C-6) 



«_ w . = _» (T+1) (i^)', 



: -3{ 



(C-7) 
(C-8) 



The solutions for F u F 3 , and are easily obtained and are 
as follows: 



Fl= U S^" {■ 4 ^" t +[('>'+l) + (3-7)/3 2 ](e- 5 «-3e- 3 «)} (C-9) 

(C-10) 



F 3 ~ 8 ^[Ase-'t+h+Wl-ffc-'t] 



F > = m [Ae_6f + (7+ 1} (1 ~ <3)<r3i] 



(C-ll) 



where A3, and .4 6 are arbitrary constants. It is noted 
that, in general, the expressions for F lt F 3 , and F s should 
each contain two arbitrary constants; however, the condition 
that F x , F 3 , and F s vanish at infinity requires that the 
omitted arbitrary constants be taken equal to zero. 

The arbitrary constants Ai, As, and A B are determined by 
means of the boundary equation (B-8). In terms of the 
variables £ and jj equation (B-8) takes the form 



(' 



sinh £ cos ij ^ — cosh £ sin 17 ^= 
(cosh £ sin 77 -^4-sinh £ cos i\ ^ (C-12) 



and the velocity potential <j> has the form 

<£=cosh £ cos y+J^ (3<r* cos 17— e~ 3 * cos 3ri) 
+<*CFi cos V +F 3 cos 3v+F s cos 5 V ) + . . . 
With 2= cos a, the boundary is given by 
y=&t sin 3 a—spp sin 3 a cos 2 a 



(C-13) 



and 



+ 9 ^^(8 sin 3 a cos 4 a— 3 sm s a cos 2 a) — 



^f= -| sin 2a+| /5<*(2 sin 2a+5 sin 4a) 



~Yg j8^(9 sin4a+14 sin 6a) + 



(C-14) 



At the boundary, inclusive of terms conta inin g the third 
power of t, 

p = l- tf+y 2 ) =sin 2 a(l - /^sin 4 a+ 6p 2 P sin* a cos 2 a) 
sin £=/3f sin 1 a(l— 3f cos 2 a) 

cosh £=1+| /S 2 ** sin 4 a-3/3 2 ? sin 4 a cos 5 a 



(C-15) 



e -f=l-/3f sin 2 a+\ sin 4 a+3/3f 2 sin 2 a cos 2 a 

sin 77= sin a ^1— | /S 2 * 2 sin 4 a+^ j^f* sin 2 a 

+3/3*^ sin 4 a cos 2 a— 3/S 2 * 3 sin 2 a cos 2 a) 
cos i7=cos a (l — I jS 2 ^ sin 4 a+SpPf 3 sin 4 a cos 2 a) 
At the boundary then, 

§ = il8 ^[-5^-3(7+D (1-iS 2 )] 

Hence, with all terms involving powers of f higher than the 
second excluded, the boundary condition (equation (C-12)) 
yields the following expressions for the arbitrary constants 
A u Az, and A*: 



(C-16) 
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4i=4[(7+D + <3— dfl-S ^Ep- 
A 3 =— g (7+l)(l— F)— g j-f^i 



4i=-|(T+l)(l-/f)+i5 



8 )S a (3+5/3+2j3 3 ) 



1-/3 2 



(0-17) 



The expressions for J\, F 3 , and then become 



^tt)=^ 1 - 1 /[(T+l)+(3-7)/3 i I(4e-«-3e^+^)-| 1 ^e- { 

3 l + fl+2,8 2 



(0-18) 



Equation (C-13), together with equations (C-18), represents a solution of the fundamental differential equation (1) and 
satisfies the boundary conditions at the surface of the solid boundary and at infinity, insofar as the terms inclusive of 
the second power in t are concerned. The present process can be extended to the higher terms in the development of 
the velocity potential fa but it can readily be seen from the complexity of the right-hand side of equation (9) that the 
labor involved would increase rapidly with the degree of approximation. In view of the importance of the problem, 
it has nevertheless been thought worth while to extend the calculation to the third approximation fa. 

APPENDIX D 

DETERMINATION OF fa 



The differential equation (9) for fa in terms of the complex 
variables s and z can be written as follows: 



1-P 
/3 s 



h(rr+ 1) (1 - (fa„+fan) (<h,*+faf) 



+ (7+I) (1-0 1 ) (*i«+*i»)*l«&I 

+ [(7+ 1) (1 - (F) +2 J 8*l[fo 1 „+&H) (fa,+fai) 

+ (*U+*ll) (&«+**»)] 

+2[( T + 1) (1 + (S 2 ) (*u+<hx)<feii 
-2/3 2 [(*u 2 -*,^) (&„-&*) + (.fa,~fai) (*i»-*ih) 



+ (.fat— <f>ii) {<hzz— fan)\ 
where, for example, 



(D-D 



&4* 
dzdz 



Again, complex variables f and f are introduced in place 
of 2 and i by means of the transformation 

2= cosh f 
According to equation (C-3), 



* l= 8^ (3e- f -c- 3 f+3e-f-e- 3 0 



(C-3) 



Consider the term 



T, =f (7+ 1) (1 - 0") -f-4je=] fan) (.faf+faf) 

+ (7+ 1) (1- P 1 ) (*i«+«ih)*u*u-2/3 s (* 1 „-*,h) (faf-faz 2 ) 



Now 



and 



, 3 g-f£ 
2/3 smhf 



4>i*+<fc?=g0i e _4f cos 4?; 
9i 

4>u—4>if= — e -4£ sin 4v 

faxfat = Jg^jl 

. 3 e~* cos 3q— e~ 3{ cos q 

* 1 " + * I "~2l sinh f sinh f 



3i e~* sin 3g— e~ 3{ sin q 
sinh f sinh f 



It follows then that 
T x sinh f sinh f= 

^ ({ [(7+ 1) (1 - ^) 2 +8^e- 6 f-2 (7+I) (1 - ^)e- 7 qcos, 
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+ {2(7+l)(l-/3 1 )e- 6£ -[(7+l)(l-/3 2 )+8i3 2 ]e- 7 «}cos 3q 
-(7+ 1) (1 - cos 5 V + (7+ 1) (1 - /S 2 ) 2 ^ cos 7tj) (D-2) 

In the handling of the terms, on the right-hand side of 
equation (D-l) that contain the derivatives of fa, it is con- 
venient to separate fa into two parts: namely, fa', which is 



(D-3) 



THE FLOW OF A COMPRESSIBLE 

a function of f plus a function of f , and fa", which is a function 
of f and X • Thus, if the variables f and f are introduced into 
the expression for fa, it follows easily that 

^=A(3e-f-e- 3 f+3e-f-e- 3 f) 

+J5(5e- 3 r-3e- e f+5e- 3 f-3e- 5 "f) 
0 a "=a4e-r-3e-f- s f4-e- SJ - 3 f+4e-f-3e- i f-?+«- 8r_lF ) 

+Z>(-8e-i-+e-r- 4 r+ef- < f-8e-f+e- 4 f-r+c- , f + f) 

where 

a 3 / i ,y i-j3* V 1-2/3 

n 9 / j-i/ 1 -^ 3 3+5/8+2^ 

D =m < T+1 >( i ^) 1 

Table II gives values of A, 5, <7, and D for values of the 
Mach number M ranging from zero to nearly unity. The 
value of y taken is 1.405. 



FLUID PAST A CURVED SURFACE 
Consider the term 

T 3 = [(y+ 1) (1 - /3 s ) +20>][(fa a +fan) (fa',+fa'i) 

+ (<£l»+ fait) (0s'»+ fa'zi)] 

— 2p 2 [(fan—<f>i3) (fa'i—fa'i) + (fat— fat) (fa' a — 4>2zl)\ 
+2[(y+l)(l + ft-20 l ](fa t +fa?)fa''* 

Now 
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and 



. 12Ar 2 r+120&;-T 
** " " inhr 

fa',+fa'l= — 12Afi-^ cos 2tj— GOiJe- 4 * cos 4j; 
fa't— fa'z=12iAe-*i sin 2i ) 4-60iBe- 4 * sin 4ij 



■feVf-fo's ^-^ ^ r [12A(<rf cos 3>j-e- 3 £ cos 
+ 120B(e~ 3 £ cos Sy—er* cos 317)] 
fa'u—fa'-5=— a - nh /sm iTf [12A(c-f sin 3ij— e" 3 « sin i?) 
+ 120B(e- 3 « sin 5i3-e" 5 * sin 3tj)] 



It follows that 

T 2 sinh f sinh T=(-j[(y+l)(l-(P)+m 46- s «+{|[( r +l)(l-/S s )+4/3 s ](24-55)+^[( 7 +l)(l+/3 8 )-2 ) 3 a ] <?)<r s * 

+(j [(7+l)(l-^)+4^]S-|[( 7 +l)(l+^)-2i3 s ](3a+2Z>)} e-^)cos v 

(T+Dd-m-f [(t+1)(1-^)+4^B+|[(t+1)(1 + ^)-2^](3C+2I>)J«- 5 « 

+{j [(7+1) (1-^+4^] 5-^[( 7 +l)(l-f-^-2^] C} e-V) cos 3, 

(7+l)(l-i3 ! )4e- 3 f+{i|^(T+l)(l-/3 J )5-^[(7+l)(l+/3 J )-2^i>J cos 5,. 

+{-^(T+l)(l-/3 ! )5+^[(T+l)(l+/S s )-2|3 a ]Z'J e -«fcos7 J , (EM) 

The calculation of Ti and 2j has not been very involved. The calculations for the remaining terms on the right-hand 
side of equation (D-l), however, are quite laborious and therefore only the final results will be presented. The terms to be 
calculated are 

T 3 =[(y+1) (l-0 r )+2$ I \(fa„+fas) (fa\+fa"i) 
Ti = — 2 0* (fau — fas) (fa" 1 — fa" 7) 
T s = [(y+1) (1- P 1 ) +2f?\ (fa,+faj) (fa" a +fa"n) 
Tt= — 20 2 (fa 1 — fai)(fa" n — fa"u) 
The expressions for these terms are as follows: 

T 3 sinh f sinh J =|[(7+1) (1— /3 s ) +2|^{2C[(2e*-3e-«-4«- < «+7«" , e-3tf-»«) cos 17 

+ (2e-«-3e- 3 «+3e- s «-e- 7 f) cos Z v ]+D[(^-B^+lQe^-Ze-^-Ae-^ cos n 
+ (4e«— 8e~«— 4e- 5 «-5e- 7 «) cos 3i}4-(4e-t— e -7 *) cos 5^4- (4te^+e~^ cos 7i?]} 

-f [(7+1) (1-^+2/S 2 ] sinh* £ cosh £ [<7(2^+«r«) +2 D]^f^j (D-5) 
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Tt sinh j- sinh f =— 6/3{2C[(2ef— 3e-*+e~ t£ ) cos Ti+.(2<r«— 3e-^+e~ 7 ^ cos 3v] 

+Z?[(4e 3 «— 8ef— 3eT^rh4e- 7 0 cos ri+(4ef— 8fl-*+4e- |S £— 5 6 - 7 «) cos 3r; 
+ (4e-«+e- 7 «) cos 5rj+ (4e- 3 «-e-*«) cos 7t,]} 

+96/5 sinh* $ cosh £ [0(2^+^+20] SS g f^ E ? (D-6) 

T 8 sinli f sinh ?=|t( T +l)(l-i^+2^{qT(-4ef+3e-«+lle- ,, «-15e- 7 f+5e- 11 «) cos T)+(-Be- 3 t+3e- s l-3e- 7 l+36- 9 e) cos 3r/] 
+2D[(— d^+4^+e-i-5e- s ^ cos u+(— e^e"- 3 *-^- 7 *) cos 3r,+e~£ cos Sr^e" 3 * cos 757]} 
(l-/3 s )+2/3 2 ] sinli 3 £ [C r (8+lle" i f+8e- 4 £+10e-«£+166- 8 £+7«- 10 «) 

+20(6^+4+2^+4^+5^1 ShT^T 

-f [(7+l)(l~i5 2 )+2i5 J ]e- s « sinh 5 £ cosh* £cosh2£ [C(2e-^+ e -^)+2l>] f (D-7) 

T 6 sinh f £anhj7=-6/3{C[(-4ee+3e-«-3e- 5 f+9e- 7 f-5tf-"^ cos rj+^e-^+Se-'S+Se-^-Se- 9 *) cos 3r;] 

+2D[(-3e 3 f+4ef-e-«+3e- 5 £) cos ij+C-ef-Se-^-Se-^cos 3r?+e-£ cos Bij+Ze^ cos 7t,]} 

-48)3 sinh<£ cosh £ [G(8e" i f+ll6-^+16e- 6 {+7e- 8 f)+2Z>(5+4e- 2 «+5e-^)] S5E C |^ r = 

+384i3 6 -^ sinh«£ cosh 3 ? [^7(2^+6-^) +2Z)] (D-8) 
Equation (D-l) can now be written as follows: 

^+^=^[We- 3 «+A 1 e- 5 «+^7 1 e- 7 «+A 1 e- 9 «+4u 1 «- u «+^- 18£ ) cos v 

+ (^ i 8 «- t «+ J i 7 8 e- 7 «+ A^+^uV^+^Ve""*) cos 3tj 

+ (^^f+^ 8 e- s «+4 7 5 e- 7 f+A 5 «^ 5 +^uV 11 «+As 6 e- ut +^i7 5 «- 17 «) cos 5r, 

+ (AV*«+4 7 7 e- 7 «+A 7 e^ f +A l 7 e- u «+^i S 7 e- w «+^i5 7 e- w «+^ w 7 e-" { ) cos 7ij] 

+^j^(f (7+1) (l-/3 s )[C(-l2e- e «+21e- 7 f-96- u e+20e- 13 «-33c- I6 £+13e- w £) 

+2 J D(-5e- s E+12e^-7e- 7 «+9e- u «-20e- ls «+lle- u f)l 

+12i3[(7(-16e^+27e- 1 «+20e- u «-44e- w f+13e- u ^+22?(-7e- 7 «+16e^-20e- 18 f+lle-"0]J 2J e_ifl£ cos (2«+7)t; 
+^^{|(7+l)(l-/3 , )[tf(-2e- 5 f+3e- 7 «-6-^ 

+24i3[C?(-2e- 9 f +3e- u «+2e- 13 «-46-"«+e- 1 ^ +2Z>(-<r 7 £ +2 e - 9 £-2 e - 13 £ +e-^]] S tt* -2 "* cos (2n+7)r; (D-9) 

J n-l 

where expressions for the constants A Q V are given at the end of this appendix. 
The right-hand side of equation (D-9) suggests a solution of the form 

<*>3=<7i(£).cos tj+tfstt) cos 3r/+<? 6 (£) cos 5t;+<7 7 (£) cos 7ij+S6 r 2»+7 cos (2»+7)r; (D-10) 

n-l 
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© 

When this expression for <fe is inserted into equation (D-9) and the coefficients of cos 17, cos 3t), . . . on both sides of the 
equation are equated, the following differential equations for Q u G s> Q s , G 7 , and G^ are obtained: - 

If 1 - <?i -^n (A^i+A^e-^+A^+A^+An'e-^+A^-^ (D-ll) 

^| 3 - 9ft =±=f (Afe-v+ASe-V+^-'t+An'e-^+Au^ (D-12) 

^ s -25&=^^ (^V 3£ +^r*+A s «" 7£ +A^ (D-13) 

^-49(? 7 -^ (AV^+^V^+A^+^LnV^f+^sV^f+^/e-^f+^/e-^ (B-14) 

^p*?- (27i+7) s (? 2a+7 =^=^j| ( 7 +l) (l-^[<7(-126- t «+21e- 7 «-9«-»f+20e- 13 «-33 e - lt f+13e- 1 »f) 

+2Z>(-5e~ 3£ + 12<r 6£ -7«- 7£ 4-9<T 11£ -20e- IS£ + 1 l«s- 15£ )] 
+ 12/3[(7(- 16e: t «+27e- 11 f+20e- 13 f-44e- 15 f + 13e- lfl f) 
+2Z> (- 7e~ 7£ + 16e- fl «-20e- 13 i+ 1 le- lt£ )])e- i,,£ 

(7+1) (l-/S I )[C , (-2e- t £+3e- 7 f-e- n £+2e- ls £-3e-^+e- lt| f) 

+2Z>(-«^e+2<r*«-e- 7 «+<5-"«-2e- 13 t+<r ls f)] 
+24/3[0 r (-2e- 9 «+3e- 11 «+2e- 13 «-4e- w «+e- 19 «) 

+2r>(-e- 7 f+2e- s «-2e- w «+e- ls ^]Jne- s »f . (D-15) 

The solutions of these equations are easily found and are as follows: 

Gi®=>^r- (ci«~ £ +|j ^ 1(! " 3{ +^ le " 5f +^ A, 1 e- 7 i+^ Afe-'t+r^ Axtf-w+JL A^e- 13 *) ■ (D-16) 

<%(£) (<^«" 3? +B ^ 3e_t£ +4^ ^ 3 «" 7£ +4 ^ 8e_S{ +n2 4 » 8e ~ llt+ 2l6 4 » 3<r " £ ) CD-17) 

GM)=-=f (c^-^Afe-*^ Ater*+±W*+± A ^+^A U ^+^A U ^+^A„^) ' (D-18) 

W (^" 7{ -^ ^ 7e ~ 5J -n ■ A77fe_7? +3^ - 4 » 7 «"' t +^ ^" 7e_U£ +TM ^ 137e ~ I3t+ 176 4 * V "" ,+ 8l2 ^ 7e ~ U£ ) C^Htf 

<W*)=^( ^ + 7+| (7+1) (1-/S*) ^-I^+^^-h*" 8 ^ 6 " 11 *) 

+2P (l^-J^+l t+Jart-J^+i e -s f )] 

+12,3 [tf (-| <T S£ +§ e-^+i e"^)+2P (| £+± ar*-l e -<*+-L e -^J j e - 0 , +nf (D _ 20) 

The expression for the velocity potential 4 can now be written as follows: 
<£=cosh £ cos v+^p (3e~ £ cos 77— e~ S£ cos Zv)+^(Fi cos y+F 3 cos 317+^ cos 5rj)+f cos i)+G 3 cos 377+^ cos 5t] 

+ G r cos 7^+^) (cosh^Z-cos"^" e ~ t * cos ' 7_e_3f 003 3 ^- c_5£ cos 5>?-«- 7£ cos 777)] (D-21) 
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where F u F 3 , and F t are given by equations (C-18) and 6 0 is the following function of £: 

+ 12 ^ [ci-^+l^+l^-Z^ (D-22) 

In this expression for 6 0 the result has been anticipated that the arbitrary constants <^ B + 7 are independent of n and may bo 
taken to be C 0 . The constants C 0 , d, Cz, O s , and d are determined by the boundary condition (equation (C-12)) 



^sinh J cos 17 -^—cosh £sinr;-|^ ^==/3 2 ^cosh £ sin 17 -^+sinh f cos rj 



(C-12) 



If the expression for <f> given by equation (D-21) is substituted into the boundary condition and if equations (C-14) and 
(0^15) are used, it follows that 

°° = ~lh (y+1) (1-j82) WC-Z4D)+\ mC-MD) 

Hb i=? [f ^+31^+23^+^+9- f (7+D (^T)V+| (7+D (*f£)V] 

-i (-f 6 ^ 5 +^ ^ 5 +<S ^+1 ^ 5 +m V) . 

_if_A lir-L . 11 A 7 7 , 2L A 7 .19. A 7^ 

7\, 24 ^ + 14^ 7 + 32^ 9 + 72^ u + 120 + 176 + 312 Al V 



The expressions for the constants A t v are as follows: 
^=§ Ky+lHl-F)+m(-ZA-4C+8n) 

+| (7+1) (1-/8*) (84-153+280-70) 



+12/3 [(37+22)0+64-155-80] 
27 



32, 



( T +l)(l-/3») 



+| (7+l)(l-/S 2 )(305-19C-24Z?) 
— 12/3 [ 7 (30+20)— 3QB+70+2O0] 

A 1 7(7+1) (1-/3*) (30-20) -120,3(20-0) 



^n 1 =7[(7+l)(l-/3 a )+4/S s ]a 

27 Q 
A 3 =g^ (7+l)(l-/S 4 )+| (7+l)(l-/3 a )(2^1-10B+5(7) 

+12/3 [ 7 (30+20) -305-50+ 120] 

^7 3 =-^ [(7+l)(l-/3 s ) i +8/3 J ] 

+| (7+l)(l-/S 2 )(15B-16(7-I>) 
+12/3 [(-37+10)0+155] 



A, s =|[(7+l)(l-|8 ! )+4/S ! ] (50-80) 



21 



A u s = -j ( T + 1) (1- /3 s ) (30-20) -168/8(20-0) 
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97 

a 16 *=j [(7+1) a-m^ g 

^3 B =-| (y+l)(l-if)(M+40-8Z7) 

A 6 =f (7+1) (l-fiC-^ [(y+l)<H-^)+m D 

97 1 e; 

A ' = ~m ( V+l)d-/3 2 ) 2 +j(7+l)(l-/S 2 )(9B-Z?) 



-12/3 [1*2(7+ (2 7 -23)2?] 
21 



~ (y+1) (1-/3*) (0-2Z?) +252/80 

Au 6 -=-| [(7+l)(l-/3 a )+4/3 s ] (0-22?) 
^.3 e =~j (7+l)(l-/5 2 )(3(7-22?)-216/S(2<7-2?) 
^n J =f [(7+1X1-^+4/3*] a 



APPENDIX E 



r 

97 

A 7 =^g3(7+l)(l-/3T 

-| (7+1) (1-/3 2 ) (455+ 12(7-292?) + 12/3(2 7 +l)Z> 

^7 7 =f (7+D(l-^)C-f [(7+D(l-/3 2 )+4/3 1 ]2? 
A 7 =-192/3(C-22?) 

-4ii 7 =-f[(7+D(l-)3 2 )- 

j [(7+l)(l-j8 J )'+4/S s ]«7-22?) 

-j [3(7+i) (l-z^+ie^c 

f [(7+D(l-|8 J )+4^]2? 
j [(Y+l)(l-ff>+4ffl < 



-12/3 2 ] (7+^ ( 7 +l)(l-^2? 



+' 
7= 39 



\G 



THE FLOW PAST A CURVED SURFACE BY THE POGGI METHOD 



Equation (D-21) represents a solution of the fundamental 
differential equation (1) that satisfies the houndary con- 
ditions at the surface of the solid body and at infinity, 
insofar as the terms inclusive of the third power of the thick- 
ness coefficient t are concerned. The method used has been 
called Ackeret's iteration process and is valid for all values 
of the Mach number M from zero to unity. On the other 
hand, the method of Poggi yields the components of the 
fluid velocity in the form of power series in M 3 . Since both 
Ackeret's and Poggi's methods provide solutions of equa- 
tion (1), the two_ solutions must agree in the region common 
to both. The flow past the shape treated in the present 
paper will be calculated by means of Poggi's method and 
compared with that obtained by Ackeret's method. 

Poggi's method consists in regarding a compressible fluid 
as an incompressible fluid with a continuous distribution of 
sources in the region external to the solid boundary. In 
order to express the intensity of the source distribution, 
it is first necessary to determine the incompressible flow, 
which serves as the zero approximation to the flow of the 
compressible fluid. The first-order effect of compressibility 
on the velocity of the fluid is then given by a set of double 
integrals extended over the entire region of flow. In 
reference 5 the surface integrals are replaced by line integrals, 
which are evaluated by the methods of the calculus of 
residues. For the example treated herein, the general re- 
sults given in reference 5 are immediately applicable and 
are as follows: 

Let Z=Z(Z') be the conformal transformation of the 
profile in the Z-plane into the circle of radius c with center 
at the origin of the Z'-plane. Then, if Wo is the complex 



velocity of an incompressible fluid past the circle, the follow- 
ing expressions are formed: 



(E-l) 



The complex velocity W of the compressible fluid past the 
profile in the Z-plane is given by 



W=(W 0 +W 1 H*)^ 



(E-2) 



where only the terms involving the square of the Mach 
number M are shown. The expression for Wi is as follows: 

w ^-M? {Z,) -z*Kt)-w>K^ z ')~\ 

--^j. s F(Z',Z')+\ U(l-^) (E-3) 

where ~F (jg;> Z'^ is obtained from f(Z',Z') by replacing 
Z' by 2^ S{Z') denotes the sum of the residues of 

2,'—Z' r ^ (^''Z*) a ^ ^ e P°^ es witbin the circle G\ of radius 
c, the subscript P being dropped after the evaluation; and 
S (j^J i 8 obtained from the expression for £(Z') by replac- 
ing Z' by -gr 
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For the present example, according to equations (A-3), 
(A-4), and (A-6), 



Z=Z' 



and 



Z' + 3Z' 3 
It follows that 



(A-3) 
(A-4) 



F{Z',z')=m 



1-^75 



+2W 



and 



F (z>, ^'"^ 



(Z'-ttf) (Z'+t<*)(z'-f) (z'+^T) 



+2Z7W 



Z' 

(Z'-«f) s (Z'+i<Z) : 



2C 2 



-Z' 



(E-5) 



Consider now the contour integral 

In the first term on the right-hand side of equation (E-5) 
only the simple poles at Z'= ±id are internal to the circle <X 
The contributions of these simple poles to the residue are, 
at Z'=id, 

1 



and, at Z'=—id, 



2id(c 2 -d 2 ) Z' P -id 

P»c* 1 
~2id(c*-<F) Z'p+id 



If the two expressions are added and the subscript P is 
dropped, the result is 

C7V 1 



(E-6) 



The second term on the right-hand side of equation (E-5) 
can be rewritten as follows: 



2C7 3 



v? 



-Z' 



The contributions of the double poles at Z'=id and Z'=—id 
to the residue are, respectively, 

IPc* Z' P -2id . l iZPc'd 1_ 

2id (Z' P -id) 3+ 2 "?=^ Z'p-id 

+i^^(ZV^ 10 ^ 

and 

77V Z' P +2id 1 iTP<?d 1_ 

2id (Z'p+id) 1 2 ?=^ZVR3 

If the two expressions are added and the subscript P is 
dropped, the result is 



2U»*P E7W 1 

"(Z /s +<Z*) a 

+|PV+J) ( f^Jy bg £ff 



CE-7) 



If the expressions (E-6) and (E-7) are added, it follows that 

S(Z0 = Z7v(l+^log^(f^ (E-S) 

From equation (E-4), 

i'YZ' ZO — 773 C^' 2 c 2 )Z ,2 Z / 
k ' y Z"(Z' 3 +d 2 )(Z' a +<2 ! ) 

If §(Z') is formed from equation (E-8) and Z' is replaced by 
^7? it follows that 

Finally, if ^(Z',Z') is formed from equation (E-9) and Z' is 
replaced by c*/Z' ', it follows that 

F^Z J-OV ^+^^j^ 

+2^^q^p(2'+^log|±|) (E-ll) 
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By means of equations (E-8), (E-9), (E-10), and (E-ll) and by replacing Z'/e, Z'/c, WJU, and cP/c* by Z', Z', W u and e, 
respectively, it follows that 

W Vi , eZ /3 -l 1 , Z /a -l 

^'"^ 2 ° g l-£A(-Z' s +6) 2i " (eZ' s +l) a J + ( € +Z /a ) (eZ /2 +l) 

The complex velocity W of the compressible fluid past the profile in the Z-plane is given by equation (E-2) where 

dZ' Z'* 

JZ~ = (Z'*—l)(Z n +e) ^ from e 1 uatio11 ( A_3 )) (E-13) 
TT 0 =1— ^5 (from equation (A-6)) (E-14) 

and TTi is given by equation (E-12). 

The velocity at the surface of the profile is obtained by means of equations (E-2), (E-12), (E-13), and (E-14) with 
Z'=e« and Z'=e-«. Thus 

^- = 7T+2T^¥+? 0+ 1+2. f {«C1— *». 2.) -fi±s£ [l- 1+2 2 £ ( c -f g+ J log 1±? 



For 0=0, 



For 0=|i 



W«.=izi(i+x L(i=^~ 1+ 2 (r^J log I^-1T^ log I=7zJ) ®~ 17) 



The ve!ocity_along the X-axis external to the profile is obtained from equations (E-2), (E-12), (E-13), and (E-14) 
by taking Z'=Z'=X'. Thus 

2 ™ x ,s +A i+ 4 L u+ ; v 2 10g i- e y (X' a + e ) 2 (eX' 3 +i) s + (X' 2 +6)(6X /s +i) 1+ (x ,a + £ ) a 

26(l-6 i )X 3 (X a +l) / T/ l+c ■ _ 1 2-/IX / M] ' ™ im 

- (X' a + £ ) a (6X' a +l) a V + 2Vl X^=7jJJ 

The velocity along the F-axis external to the profile is obtained from equations (E-2), (E-12), (E-13), and (E-14) 
by taking Z'=iF' and Z'=-iF'. Thus " 

F /a l , m 2 r . . / l+e, i+e \ € r /6 +(e 8 -46+i)F /4 + e z /2 F /a , , r" 

ffr-o F 5= "el i " t "4L U " 1-e; ^ i 2 10g l- e y (6-F ,a ) a (l- £ F ,a ) 2 (e- F' 2 ) (1— £ F ,a ) ~ 1 + ( £ - F /2 ) a 

+2^ (V-ii- 6 ioe ^ a +f) (r 8 +i)'-2(i+j»rn 1 

+2er r 2V* g FWe7 (F' 4 +i)(i- e r' 2 )( e -F' i ) 4 Jt CE- 19 * 

Table III gives values of the velocity corresponding to points along the profile for the numerical case e=y (or £=0.10) 

at Mach numbers M=0.50 and 0.75. Figure 1 shows a comparison of the velocity distribution for M=0.50 with the 
corresponding calculation according to the iteration method. 

It is noted that by Poggi's method the solution of the problem is given by the components of the fluid velocity, 
whereas by Ackeret's method the solution is obtained in the form of the velocity potential of the fluid motion. Before 
a comparison of the two methods can be made, therefore, it is necessary to obtain from the velocity potential (of Ack- 
eret) the velocity components along the coordinate axes. This calculation is performed in appendix F. 
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APPENDIX F 

DETERMINATION OF THE VELOCITY COMPONENTS 



In terms of the complex variables f and f the velocity 
components in the direction of the coordinate axes are given 
by 

1 b<f> . 1 d0 



b<f> 



dX~ sinh f St R inh f df 



f 2>f sinh 



1 

inh f &]7 



(F-l) 



where the velocity potential ^ is obtained from equation 
(D-21). According to equations (6) 



»=ifli+i 2 c J +t 3 0 3 + • • • 
where, in general, 

In terms of the complex variables £ and f , 



(F-2) 



Therefore, 



3 

«i=— 2^ e -s * cos 2r/ 
3 

t>i=— 7j e -2 * sin 2 77 



(F-3) 



Similarly, with reference to equations (D-3), it follows that 

«,=4C , (2-3e^+3e^-« HI «)+8D(«*-2+«- 2 e) 

+4(4D-ZAe- 2 t+4De-*£) cos 2, 

+4[4Z)e- 2 «+(D-15B)e-^ cos 4ij 

-32 sinh 3 g cosh g[2P+g(2g-«*+g-q)] 

cosh 2£— cos 2ij [ (F-4) 

« s =4i8 sin 2, j(4I>-3^e-^+4Z>e^f) 

+2[4De- 2 t-(D+ lSB)(T*e[ cos 2ij 

4 sinh 2 tpD+CQe-^+e- 4 *)] 
cosh 2|— cos 217 

The general expressions for the velocity components U3 and 
Vs are too cumbersome to be given here. Instead, only the 
expressions for the velocity components along the profile will 
be given. Thus, along the profile, if powers of t higher than 
the third are neglected, 



u 3 =(Q l -\-%G z +5G b +7G 7 ) 0 +2{ZG 3 +5Q l ,+7Gi\ cos 2a 

+2(5G f 6 +7(? 7 )o+14(G' 7 ) 0 cos 6a 
C ! =-2 J 8t(G' 3 , +(? e , +(?/) 0 sin 2a 

+ (G/+G/) 0 sin 4a+(G/) 0 sin 6a] 



CF-5) 



where the primes denote differentiation with regard to the 
independent variable £ and the zero subscripts denote evalua- 
tion for f=0. Explicit expressions for (Gi) 0 , (#3)0, (G 5 ) 0 , (<? 7 )o, 
(C?i')o, (G 3 ') 0 , {Qt')o, and f(?/)o are given at the end of this 
appendix. 



At the boundary then, if terms containing powers of t higher than the third are neglected, equations (F-2) become 
«=l+[^ (7+D Q^J+l ?+[tq (2+p) + (G 1 +3G 3 +5G,+7G 7 ) o y 

+|-|s t+ h *+\]k 4+2 ^ +3 ^ +2(3g 8 +5g t +7g 7 ) 0 ] *» } cos 2a 
+{ W+i ^] 4 -™f+2(^ s +7. 7 )o] , } cos 4. 

32 



+1 32 [ — £ +14((? 7 ) 0 j i 3 cos 6a 



»=[-§ *+§ ^+3^ (4-2/3+5 i 3 2 )< 3 -^ (t+1) < 3 -2^(G'3 / +G r s '+G r 7 , )o< 3 ] sin 2a 

+[| <2_ A (4+6/3+3/3 i )/ 3 +^ (7+I) f 3 — 2 J 8(Gf 6 '+^/)o* 3 ]] sin 4a 

+[^ (12+260+7/3*)-^ ( 7 +l) ^^- i -2/3(G'/)o] * sin 6a 

In order to compare the velocity of the fluid given by equations (F-6) with that obtained by Poggi's method given 
by equation (E-15), explicit expressions for equations (F-6), with powers of the Mach number higher than the second 
neglected, are obtained and may be expressed as follows: 



(F-6) 
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«=(l-§ P+Jq « 3 )+(-| *+§ ?) cos 2«+(3^-g cos 4«- 2 ^ * cos 6« 

^-TM <3 )+(-t '+1 ^-|§§ ^) cos 2 «+(l '"I 0 «« ' cos 6a] 

(F-7) 

"=(-| *+f <2 +§ *") 8111 2a +( 3fS -ii <*) ^ 4a -^ * 6a 

+ [^f 3 sin2a+^f 2 -|^sin4a-|< 3 sin6 a J 
The magnitude g of the velocity, to the same degree of approximation, is then given by 

<K 1- ^ <a ) + (~l <+ ! f ) cos 2a +(i M 0 003 4a_ ^ ^ 009 6a 

+^° [(! '-wo <+ i <2 -iS ') cos 2a+ (r6 0 «» *-S 43 cos 6 «] - ^ 

which is to be compared with the corresponding expression obtained by Poggi's method. Thus, if in equation (E-15) e is 
replaced by 3f/(2+£) and all terms are expanded in powers of t up to and including f, the resulting equation is as follows: 

W<fc=(l+^ P-Jq *)-(! *~l P+^Hi <") cos 2d+jt(?-1?) cos 40-^ f cos 60 

+M » [(J f-g *)+(-§ I+| f-^g- ?) cos 20 +^(l^f) cos 40-f| f cos 6tf] (F-9) 

From the first of equations (A-9) with X=cos a, it follows that 

cos a—(l — t sin 2 0) cos 0 

and therefore 

cos 2a=cos 26— sin* 20^f— | <* sin s 0^ 

cos 4a=cos 4:0— 2t sin 40 sin 20 
cos 6a=cos 60 

If cos 2a, cos 4a, and^cos 6a in equation (F-8) are replaced by these expressions, equation CF-9) is again obtained. 
The expressions for (Qi)q, (03)0, (Os)o, (.0?)o, (,6i')o, (6 r 3 / )o, (Gs')o, and (G/) 0 are as follows: 

(^[i,-i4-^4^Hf)^-i^i)(f)l] ■ 

+ (7+1) (-§ A+f 2?-| <7+| 2?) +/s(-64+15B-~ <7+14p) ] 

»4[-f^^c«>(f)''4w)(D'g 

-f[db<*+ 1 >< 1 -«-510<*+ 1 >f-8^ 

< e ')«=iio[f ^+"+f+^+|-|h+«(f)'*+l^+ 1 )(f)*)] 

-f[-ilo^+«f-!^+^+<(-l^+!^fo«)->] 
^')"-b[ _ I is_ t-|-^~^ + I (i ' +i) (t) "~i (T+I) (^)?] 

~f [i4^+ 1 'f+f^ D -S^+ 1 »T (45B - 17D)+ ^ I> ] 
(ft0.-a [-^+1+5+^+5 ftr+D (f)'<»+| (t+» 
W ,.4[^ + 5-|-^-|(, +1 ,(f)- ( ,-|(, +1 ,(f) , i] 

w)^g,+f+§-j+i+^-i( T +i)(f)',+i(,+i)(f) , a 



<3 
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TABLE I 

VELOCITY AND PRESSURE DISTRIBUTIONS AT THE 
SURFACE OF A BUMP, «-=0.10 
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X 


Y 


Q 


c r 
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(Equation 




CE^uatlon 




CA-10)) 


(A-ll)) 


0 


LOOO 


0 


a 8760 


0.2344 


s 
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80 


0 
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TABLE II 

VALUES FOR A, B, C, AND D OBTAINED FROM 
EQUATION (D-3) 



M 
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A 


B 


C 


D 


0 


1 


0.06250 
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0 


0 


.10 


.99499 


.06221 


—.06285 


.00095 


0 
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.037a 
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.34661 
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6.81182 
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.28000 


3L0576 


L 97342 


8.8910 


11.653 
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.19900 


132.428 


9.30323 


34.428 


49.732 



TABLE m 

VELOCITY DISTRIBUTION AT THE SURFACE OF A BUMP. 
i=0.10, ACCORDING TO THE POGGI METHOD 
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X' 


y 


X 


Y 
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ff 
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M- 
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M- 
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.3118 


.0830 


1.117 


.0645 


1.133 


"1.153 






.9 


.8965 


.8676 


.8444 


SO 


.1737 


.9848 


.1668 


.0965 


1.153 


.0795 


1.173 


L198 






.976 


.8704 


.8258 


ireio 


90 


0 


1.000 


0 


.1000 


1 167 


.1002 


1.192 


1.223 






L0 


.8568 


.8072 


.7679 



TABLE IV 

VALUES OF (<7i) 0 . (GW 0 . (GWo. (<Wo GIVEN AT END OF 

' APPENDIX F 



M 


0 


(ffi) 0 


(Gi), 




(Ol), 


a so 


a 86603 


0. 10344 


-0.79011 


1. 6714 


-a 91449 


.75 


.65144 


2.3075 


-2.1109 


18976 


-1.9513 


.83 


.56776 


9.8911 


-6.7071 


6.9329 


-4.0669 


.90 


.43589 


68.849 


—12.316 


26.227 


-19.649 



TABLE V 

VALUES OF a,, a,, AND a, OBTAINED FROM EQUATIONS (14) 



\ 


M 






I 








0.60 


- 0.75 




0.90 


X 


\ 






| 0.83 


1 


ai 




0 


1.7321 


2.3028 


2.6893 


3.4412 




.1 . 


L6974 


2.2666 


2.6356 


3.3724 




.2 


1.6936 


2.1184 


2.4742 


3.1659 




.3 


1.4203 


1.8881 


2.2052 


2.8218 




.4 


L1778 


1.6658 


1.8287 


2.3400 




.6 


.86603 


L1613 


1.3447 


1.7205 




.6 


.4S497 


.64473 


.76301 


.00365 




.7 


.03464 


.04606 


.06379 


.06882 




.8 


-.48497 


-.64473 


-.76301 


-.96366 




.9 


-1.0739 


-1.4276 


-L6674- 


-2.1338 




.975 


-L6610 


-2.0762 


-2.4237 


-3. 1014 




LO 


-L7321 


-2.3026 


-2.6893 


-3.4413 


ai 




0 


2.2743 


6.4350 


10.138 


19.263 




.1 


2.0399 


4.9892 


9.3921 


18.070 




.2 


1.3673 


3.7088 


' 7.2538 


14.618 




.3 


.34809 


1.7641 


4.0016 


9.4302 




.4 


-.86531 


-.66063 


.10228 


3.1062 




.6 


-ZD593 


_-2.8678 


-3.7000 


-3.1035 




.6 


-2.9594 


-4.6469 


-6.8342 


-8.0406 




.7 


-3.2326 


-6.2746 


-8.0106 


-ia022 




.8 


-2.4744 


-4.0011 


-6.0811 


—7.0773 




.9 


—.23520 


.01270 


.33212 


2.9801 




.978 


2,7276 


6.4013 


8.0381 


10.723 




L0 


4.0063 


7.7376 


12.825 


22.703 


a> 




0 


3.4379 


21.877 


72.397 


437.84 




.1 


2.6124 


18.607 


63.333 


386.94 




.2 


.07542 


9.4933 


38.871 


261.93 




.3 


-2.9313 


-2.1171 


6.6216 


72.788 




.4 


-6.1915 


-12.039 


-22.646 


-92.633 




.5 


-6.8962 


-16.837 


-37.884 


-184.07 




.6 


—2.7710 


-11.097 


-31.791 


-104.69 




.7 


2.2311 


1.4468 


-6.9844 


-41.666 




.8 


7.0916 


14.258 


22.616 


92.680 




.9 


1 2652 


10.459 


14.608 


23.200 




.975 


-.82172 


-2a 694 


-66.978 


-386.70 




LO 


-10.033 


-39.881 


-101.38 


.-633.92 



TABLE VI 

VELOCITY DISTRIBUTION FOR A BUMP, <=0.10, CAL- 
CULATED BY MEANS OF TABLE V AND EQUATION (13) 



THE FLOW OP A COMPRESSIBLE 



FLUID PAST A CURVED SURFACE 



327 



TABLE VII 

VALUES OF CRITICAL VELOCITY OBTAINED FROM EQUA- 
TION (18) AND MAXIMUM VELOCITIES FOR A BUMP, 
*=0.10, OBTAINED FROM EQUATIONS (17) AND (E-17) 





g. 


ma 




M 


















method 


method 




0 


1.167 


L167 


CO 


.2 


1.171 


L171 


4.678 


.3 


1.176 


L176 


3.067 


.4 


1.185 


1.183 


2.316 


.6 


1.199 


1.192 


L869 


.6 


1.222 


1.203 


1.674 


.7 


1.286 


L216 


L366 


.8 


1.371 


1.231 


1.212 


.86 


1.623 


1.239 


L145 


.90 






1.093 


.98 






1.044 


1.0 






1.000 









TABLE VIII 

VALUES OF qlq„ FOR A BUMP, <=0.10, WITH Af =0. 83 



X 


7 


tiler 


0 


L443 


1 230 


.1 


1.421 


1.211 


.2 


L369 


1.160 


.3 


L267 


1.080 


.4 


1.169 


.9881 


.6 


L069 


.9028 


.6 


.9762 


.8314 


.7 


.9193 


.7837 


.8 


.8866 


.7668 


.9 


.8444 


.7200 


.976 


.7910 


.0744 


L0 


.7679 


.6481 



TABLE IX 

COMPARISON OF THE PRESSURE DISTRIBUTION AT THE 
SURFACE OF A BUMP, i=0.10, FOR Af=0.83 OBTAINED 
BY MEANS OF THE ITERATION.THE PRANDTL-GLAUERT, 
AND THE VON KARMAN METHODS 



X 






Iteration 
method 
(Equation (20)) 


Prandtl- 
Qlanert 
method 


von Karm&n 
method 


0 


-a 9133 


-a 6342 


-0.7376 


-a 3637 


.1 


-.8877 


-.6142 


-.7107 


—.3426 


.2 


-.7389 


-.6668 


-.6337 


-.3100 


.3 


-.6188 


-.4666 


-.6202 


-.2602 


.4 


-.3294 


-.3614 


-.3810 


-.1960 


.5 


-. 1184 


-.2160 


-.2267 


-.1199 


.6 


.0614 


-.0647 


-.0666 


-.0361 


.7 


.1617 


.0786 


.0762 


.0426 


.8 


.2221 


.2062 


.1972 


.1160 


.9 


.2868 


.3196 


.2986 


.1783 


.976 


.4048 


.3936 


.3620 


.2196 


LO 


.4708 


.4163 


.3812 


.2322 



TABLE X • 

VALUES OF a,, a,, AND a, CALCULATED FROM EQUATION (17) 



M 


«i 


ai 


a 


0 


2.60000 


1.60000 


1.50000 


.2 


1.63093 


1. 68726 


L67024 


.3 


L67243 


1. 71149 


1.93367 


.4 


1.63663 


L 92091 


2.43349 


.6 


1.73205 


2.27426 


3.49200 


.6 


1.87600 


2.91036 


6.76400 


.7 


2.10042 


4.22402 


12.66611 


.8 


2.60000 


7.76704 


43.7829 


.86 


2.84747 


12.52960 


112.7760 


.90 


3.44124 


26.7418 


446847 


.92 


3.82733 


38.8963 


964.722 



TABLE XI 

MAXIMUM VALUES OF THE PRESSURE COEFFICIENT 
C„u CALCULATED BY MEANS OF EQUATION (20) 



-(Cr.ii ) ... 


N, I 


0.05 


0.05 


0.08 


0.10 


0.12 


0.12 


0.15 


M \ 




(a) 












0 


0. 16406 


a 16406 


a 27744 


a 36000 


a 44856 


a 44856 


a 69344 


.2 


.16764 


. 16773 


.28369 


.36330 


.45913 


.45994 


. 60787 


.3 


.17248 


.17266 


.29221 


.37964 


.47363 


.47537 


. 62779 


.4 


.18008 


. 18035 


.30570 


.39771 


.49688 


.49963 


.66002 


.6 


.19166 


.19188 


.32661 


.42606 


.63376 


.63857 


.71196 


.6 


.20971 


. 20937 


.36000 


.47200 


.69437 


.59401 


.79907 


.7 


. 24055 


. 23753 


.41994 


.55715 


.71014 


. 69006 


.97236 


.8 , 


30606 


.28925 


.66134 


. 77066 


1.0165 


.87903 


L4634 


.86 


! 38190 


! 33621 


! 74888 


1.0737 


L4759 


L0664 


22485 


.90 


.59091 


.42108 


L3551 


2.1252 








.92 


.82480 


.47966 














0.18 


0.18 


0.20 


0.21 


0.22 


0.25 


0^5 


\ 




(a) 










(a) 


0 


a 75384 


a 76384 


a 87000 


0.93098 


a 09396 


1.1953 


1.1953 


.2 


.77275 


.77541 


.89226 


.96504 


L0199 


1.2274 


1.2352 


.3 


.79899 


.80489 


.92326 


.98859 


1.0561 


1.2724 


L2803 


.4 


.84179 


.85192 


.97408 


1.0437 


1.1158 


1.3470 


1.3793 


.6 


.91186 


.92436 


1.0581 


1.1363 


1.2164 


1.4732 


L5208 


.8 


L0316 


1.0403 


L2031 


1.2943 


1.3891 


1.6S84 


1.7666 


.7 


1.2782 


1.2431 


1.6088 


1.6326 


1.7623 


2.1889 


2 1973 


.8 


2.0176 


1.6781 


2.4539 


2.6838 


2.9491 


3.8118 


3.3117 



» M ethod of Ton Kannan. 



TABLE Xn 

CRITICAL AND LIMITING VALUES OF M AND CRITICAL 
VALUES OF C P ,ir CALCULATED BY MEANS OF 
EQUATION (21) 



M 


-(C„„) „ 


( 


-Co 








a 45 


2. 76639 


0.05 


a 16406 


0.832 


a 880 


a so 


.50 


2.12953 


.08 


.27744 


.776 


.865 


.77 


.65 


L 66519 


.10 


.36000 


.742 


.833 


.92 


.60 


1.29100 


.12 


.44856 


.712 


.816 


L09 


.66 


1.00661 


.16 


.69344 


.670 


.790 


1.33 


.70 


.77768 


.18 


.75384 


.834 


.760 


1.68 


.75 


.69008 


.20 


.87000 


.610 


.743 


1.74 


.80 


.43381 


.21 


.93098 


.698 


.736 


1.82 


.85 


.30124 


.22 


.99396 


.687 


.725 


1.90 


.90 


.18605 


.26 


L1963 


.658 


.698 


2.16 


.95 


.08783 












1.00 


0 













TABLE XIII 

VALUES OF THE PRESSURE COEFFICIENT (C.j,)^. 
CALCULATED BY MEANS OF EQUATION (22) 



Af 


-{.C„id.u 


M 




a 70 


2.90508 


L25 


a 91103 


.76 


2.53064 


L30 


.84230 


.80 


2.22420 


L35 


.78106 


.86 


L97022 


L40 


.72627 


.90 


1.75739 


3.46 


.67705 


.96 


L67727 


1.60 


.63266 


LOO 


L 42349 


1.66 


.69260 


LOS 


L 29116 


1.60 


.55606 


LIO 


L 17844 


1.65 


.62286 


L15 


L 07636 


L70 


.49266 


1.20 


.98863 


1.76 


.48481 



